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5 ROBUST PREDICTIVE DECONVOLUTION SYSTEM AND METHOD 

Technical Field 

This invention relates to a method and system for predictive deconvolution, which is otherwise 
10 known as pulse compression in radar applications. More particularly, the invention relates to robust 
predictive deconvolution using minimum mean-square error reiteration. 

Background Art 

15 In many sensing applications, such as radar pulse compression, sonar, ultrasonic non-destructive 

evaluation for structural integrity, biomedical imaging, and seismic estimation, it is desirable to estimate the 
impulse response of an unknown system by driving the system with a known signal having a finite temporal 
extent. The process of separating the known signal from the received output of the system in order to 
estimate the unknown impulse response is known as predictive deconvolution. 

20 Predictive deconvolution provides a means to obtain the high spatial resolution of a short, high 

bandwidth pulse without the need for very high peak transmit power, which may not be feasible. This is 
accomplished by transmitting a longer pulse that is phase or frequency modulated to generate a wideband 
waveform. The transmission of the wideband waveform into the unknown system results in a received 
return signal at the sensor that is the convolution of the waveform and the system impulse response, which 

25 possesses large coefficient values at sample delays corresponding to the round-trip travel time of the 
transmitted waveform from the sensor to a significant reflecting object (also called a scatterer) and back to 
the sensor. The purpose of predictive deconvolution is to accurately estimate the unknown system impulse 
response from the received return signal based upon the known transmitted waveform. 

A well-known approach to predictive deconvolution, used extensively in radar and biomedical 

30 imaging applications, is known as matched filtering, e.g. as described in M.L Skolnik, Introduction to Radar 
Systems , McGraw-Hill, New York, 1980, pp. 420-434; T.X. Misaridis, K. Gammelmark, CH. Jorgensen, N. 
Lindberg, A.H. Thomsen, M.H. Pedersen, and J.A. Jensen, "Potential of coded excitation in medical 
ultrasound imaging," Ultrasonics, Vol. 38, pp. 183-189, 2000. Matched filtering has been shown to 
maximize the received signal-to- noise ratio (SNR) in the presence of white Gaussian noise by convolving 

35 the transmitted signal with the received radar return signal. One can represent matched filtering in the 
digital domain as the filtering operation 

40 where x MF (?) , for £ = 0, • • - , L - 1 , is the estimate of the delayed sample of the system impulse response, 
s = Ui s 2 - ■■ s N ] T is the length-N sampled version of the transmitted waveform, 
y(?) = [y(£) >(^ + l) •** y(£ + N-l)] T is a vector of contiguous samples of the received return signal, and 
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5 (•)" and (•Y are the conjugate transpose (or Hermitian) and transpose operations, respectively. Each 
individual sample of the return signal can be expressed as 

vO?) = x r 0?)s + vO?), (2) 
where x(i) = [x(£) x(E-i) x(£ - N + l)] r consists of samples of the true system impulse response and 
10 v(£) is additive noise. The matched filter output can therefore be written as 



x MF (e) = s H A T {£)s+s H v(t), 



(3) 



15 



where v(*) = [v(£) v(£ + 1) - - • v(£ + N - l)] r and 



A(0 = 



x(t) x(M) 
x(£-l) x(£) 



x{e+N-iy 



x(e-N+i) 



x(e-\) x(e) 



(4) 



is a collection of sample- shifted snapshots (in the columns) of the impulse response. 

From (4), it is obvious that estimation via matched filtering will suffer from spatial ambiguities 

20 (also known as range sidelobes in the radar vernacular) due to the influence from neighboring impulse 
response coefficients (i.e. the off-diagonal elements of A(£) ). To alleviate this effect, Least-Squares (LS) 
solutions have been proposed, e.g. in M.H. Ackroyd and F. Ghani, "Optimum mismatched filters for 
sidelobe suppression," IEEE Trans. Aerospace and Electronic Systems, Vol. AES-9, pp. 214-218, March 
1973; T. Felhauer, "Digital Signal Processing for Optimum Wideband Channel Estimation in the Presence of 

25 Noise," IEE Proceedings-F, Vol. 140, No. 3, pp. 179-186, June 1993; S.M. Song, W.M. Kim, D. Park, and 
Y. Kim, "Estimation theoretic approach for radar pulse compression processing and its optimal codes," 
Electronic Letters, Vol. 36, No. 3, pp. 250-252, February 2000; B. Zrnic, A. Zejak, A. Petrovic, and I. Simic, 
"Range sidelobe suppression for pulse compression radars utilizing modified RLS algorithm," Proc. IEEE 
Int. Symp. Spread Spectrum Techniques and Applications, Vol. 3, pp. 1008-1011, September 1998; and T.K. 

30 Sarkar and R.D. Brown, "An ultra-low sidelobe pulse compression technique for high performance radar 
systems," in Proc. IEEE National Radar Conf., pp. 111-114, May 1997. LS solutions decouple neighboring 
impulse response coefficients which have been smeared together due to the temporal (and hence spatial) 
extent of the transmitted waveform. The LS solution models the length-(L+/V-l) received return signal as 



35 



y =Sx+ v, 



(5) 



where x = [x(0) x(l) -» x(L-l)] r 

are the L true impulse response coefficients that fall within the data 
window, v = [v(0) v(l) ■-- v(L+ N - 2)] T are additive noise samples, and the convolution of the transmitted 
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waveform with the system impulse response is approximated as the matrix multiplication 



Sx = 



s { 0 ••• 0 

: s x 

s N I ••. 0 

0 s N s } 

0 0 s N 



(6) 



The LS model of (6) is employed extensively in radar pulse compression, seismic estimation, e.g. as 
described in R. Yarlaggadda, J.B. Bednar, and T.L. Watt, "Fast algorithm for € p deconvolution," IEEE 
Trans. Acoustics, Speech, and Signal Processing, Vol. ASSP-33, No. 1, pp. 174-182, Feb. 1985; and 
ultrasonic non-destructive evaluation, e.g. as in M.S. O'Brien, A.N. Sinclair, and S.M. Kramer, "High 
resolution deconvolution using least-absolute-values minimization," Proc. Ultrasonics Symposium, pp. 
1151-1156, Dec. 1990; and D.-M. Suh, W.-W. Kim, and J.-G. Chung, "Ultrasonic inspection of studs (bolts) 
using dynamic predictive deconvolution and wave shaping," IEEE Trans. Ultrasonics, Ferroelectrics, and 
Frequency Control, Vol. 46, No. 2, pp. 457-463, Mar. 1999. The general form of the LS solution is 



(7) 



For the received signal model of (5), it can be shown that the LS solution of (7) is optimal in the 
mean-square error (MSE) sense when the additive noise is white. However, upon further inspection one 
finds that the LS received signal model does not completely characterize the received return signal because it 
does not account for the convolution of the transmitted waveform with impulse response coefficients x(£) 
prior to ^=0. The result is that the presence of a significant impulse response coefficient within AM 
samples prior to x(0) can cause severe mis-estimation of the desired coefficients within the data window. 

There is, therefore, a need for a predictive deconvolution system with improved robustness, 
accuracy, and resolution. 

Disclosure of the Invention 



According to the invention, a method for processing a received, modulated pulse that requires 
predictive deconvolution to resolve a scatterer from noise and other scatterers includes (a) receiving a return 
signal; (b) obtaining L + (2A/-l)(N-l) samples y of the return signal, where y(£) = x T (£) s + v(£) ; (c) 
applying RMMSE estimation to each set of N contiguous samples to obtain initial impulse response 
estimates [x,{-(Af -l)(AM )}, ••, £,{-1}, £,{()}, ••,x ! {L-1},jc i {L},- -,x 1 {L-1-KM-1X^/-1)}] ; (d) computing power 
estimates p x (£) = \x l (e)\ 2 for £ = Af-1), - , L-l-KM-l)(Af-l) ; (e) computing MMSE filters 

according to w(£) = p(i) (C(£) + R) -1 s , where p(£) = \x(£)\ 2 is the power of x(E) , and 
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5 R = E[\(i) \ H (£)] is the noise covariance matrix; (f) applying the MMSE filters to y to obtain 

[Jt 2 HA*-2)(N-l)V^i^ ; and repeating (d)-(f) for 

subsequent reiterative stages until a desired length-L range window is reached, thereby resolving the 
scatterer from noise and other scatterers. 

Also according to the invention, a radar receiver system includes a receiver, a processor including a 
10 RMMSE radar pulse compression algorithm, and a target detector. 

Also according to the invention, a method for processing a received, modulated radar pulse to 
resolve a radar target from noise or other targets includes (a) receiving a radar return signal; (b) 

obtaining L + (2M-1)(AM) samples y of the radar return signal, where y(£) = x T (i) s + v(^) ; (c) 
applying RMMSE pulse compression to each set of N contiguous samples to obtain initial radar impulse 
1 5 response estimates [*,{-( Af-1)(AT-1)}, • - • , jc,{-1}, ^{0}, - - • , ^{L-l}, jc,{L}, - - , Xl {L-l+{M -l)(AHL)}] ; (d) 

computing power estimates p x (t) = |Jc,(€)| 2 for t = -(M-l)( • • • , L-1+(M-1)( Af-1) ; (e) computing 
range-dependent filters according to w(£) = p(£) (C(^)+ R ) _1 s , where p(£) = \x(£)\ 2 is the power of x(£) , 
and R = E[y(t) v H (0] is the noise covariance matrix; (f) applying the range-dependent filters to y to 
obtain [£J-(Af-2)(AM)},.. s ^ ; and repeating (d)-(f) 

20 for subsequent reiterative stages until a desired length-L range window is reached, thereby resolving the 
radar target from noise or other targets. 

The Reiterative Minimum Mean-Square Error (RMMSE) predictive deconvolution system and 
method of the invention provide high-fidelity impulse response estimation compared to other systems such 
as matched filtering and LS. In one embodiment, the RMMSE estimator reiteratively estimates the MMSE 

25 filter for each specific impulse response coefficient by mitigating the interference from neighboring 

coefficients that is a result of the temporal (i.e. spatial) extent of the transmitted waveform. The result is a 
robust estimator that adaptively eliminates the spatial ambiguities that occur when a fixed receiver filter is 
used. 

The invention has obvious applications in military radar and sonar, and it is also useful for civilian 
30 radar, e.g. in airport radar systems and in weather radar systems, e.g. Doppler radar, and other environmental 
radar applications. It may also find use in range profiling, image recognition for Synthetic Aperture Radar 
(SAR) and Inverse SAR (ISAR), remote sensing, ultrasonic non-destructive evaluation for structural 
integrity, seismic estimation, biomedical imaging, inverse filtering of optical images, or any other 
application requiring robust deconvolution of a known waveform (or filter) from a desired unknown impulse 
35 response. 

Additional features and advantages of the present invention will be set forth in, or be apparent from, 
the detailed description of preferred embodiments which follows. 

Brief Description f the Drawings 

40 



FIG. 1 is a schematic diagram of a predictive deconvolution system according to the invention. 
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5 

FIG. 2 is a block diagram of a 3-stage RMMSE predictive deconvolution algorithm according to the 
invention. 

FIG. 3 is a graph showing an impulse response with a single large coefficient (scatterer) in the data window 
10 according to the invention. 

FIG. 4 is a graph comparing the performance of the RMMSE algorithm according to the invention with LS 
and matched filter systems for a single-scatterer impulse response. 

15 FIG. 5 is a graph as in FIG. 3 but showing an impulse response with a second large coefficient just prior to 
the data window according to the invention. 

FIG. 6 is a graph comparing the RMMSE, LS, and matched filter systems for a scatterer just prior to the data 
window. 

20 

FIG. 7 is a graph comparing the RMMSE, LS, and the matched filter systems for a realistic scenario in 
which there are several scatterers with substantially varying magnitudes. 

25 

Best Mode for Carrying Out the Invention 

Definitions: The term "convolution" means the process that yields the output response of an input to a 
linear time-invariant system, such as is described and defined in J.G. Proakis and D.G. Manolakis, Digital 

30 Signal Processing:Principles, Algorithms, and Applications , 3rd Ed., pp. 75-82, Prentice Hall: Upper Saddle 
River, NJ (1996), incorporated herein by reference. The term "deconvolution" as used herein means the 
process that given the output of a system determines an unknown input signal to the system. See Id. at p. 
355, incorporated herein by reference. The term "scatterer" means something in the path of a transmitted 
waveform that causes a significant reflection (relative to the noise) back to the receiver of the sensor. 

35 Referring now to Figure 1, a predictive deconvolution system 10 includes a transmitter 12 for 

transmitting a phase or frequency modulated pulse (or waveform) 14 that, upon interacting with its 
transmission environment's unknown impulse response 16 (to be estimated, as is described further below) 
and forms a signal 18 that is a convolution of waveform 14 and impulse response 16, an analog front-end 
receiver 18 for receiving signal 16, an optional analog-to-digital (A/D) converter 20, a processor 22 for 

40 processing analog or digital signal 16, and a detector 24. The processor 22 includes a Reiterative Minimum 
Mean-Square Error (RMMSE) estimation algorithm, as is described below. 

Minimum Mean-Square Error (MMSE) estimation is a Bayesian estimation approach that employs prior 
information in order to improve estimation accuracy. The exact form that the prior information will take 
will be addressed shortly. First, however, the signal model must be constructed. From (3), we see that the 
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5 collection of N samples of the received return signal can be expressed as 

y{e) = A T (e)s+y(e). (8) 

This is the received signal model used by the matched filter formulation and takes into account all the 
10 necessary impulse response coefficients for a given received return sample. To develop the MMSE filter, 
the matched filter s H in (3) is replaced with the MMSE filter, denoted vt H (£), in which the form of the 

MMSE filter is dependent upon the length-TV sample swath of the impulse response under consideration. 
Thereafter, the MMSE cost function 

15 V(0=£[|x(0-w w (£)y(o| 2 ] (9) 

is solved for each impulse response coefficient € = 0,--,L-l, where £[•] denotes expectation. This is 
done by differentiating with respect to w w (i) and then setting the result to zero. The MMSE filter is found 
to take the form 

20 w(0 = (4y(0 y"(0]) _1 £[y(0 **(*)]. (10) 



25 



where (•)* is the complex conjugate operation. After substituting for y(0 from (8) and assuming that the 
impulse response coefficients are, in general, uncorrected with one another and are also uncorrected with 
the noise, one obtains 



(11) 



30 



where p(£) = E[ [x(£)| 2 ] is the expected power of x(£) , and R = E[y(£) v w (£)] is the noise covariance 
matrix. On assuming neighboring coefficients are uncorrected, the (/, j) th element of the matrix C(£) is 



*M7 

Y,P( e -» + <-!) s(n) s\n - i + j) 



(12) 



in which K L = max{0, i - j} is the summation lower bound and K v = mm{N - 1, N - 1 + 1 - j} is the upper 

bound. Also, any prior information regarding the noise can be employed via the noise covariance matrix R . 

35 For instance, for a white noise assumption R is diagonal. 

In its current state the MMSE filter is a function of the powers of the surrounding impulse response 
coefficients, which in practice are unavailable. This lack of prior knowledge can be taken into account by 
setting all the initial coefficient estimates equal. Therefore, the initial MMSE filter reduces to the form 
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w=C"'s 



(13) 



where the noise term is assumed negligible and 



10 



X s(n)s*(n-i + j) 



n = fC L 



(14) 



is invariant to the sample delay t . The initial MMSE filter can therefore be pre-computed and then 

implemented in the same way as the traditional matched filter. The inclusion of the matrix C serves to 
provide a "local" LS initial estimate. 
15 Figure 2 illustrates the Reiterative MMSE (RMMSE) predictive deconvolution algorithm for three 

iterations. In general, the RMMSE algorithm operates as follows: 

1) Collect the L + (2M-l)(N-l) samples of the received return signal 
Ly{-(M-l)(AM)} r -^;y{-l},^ which comprise the length-L 

data window along with the (M-l)(yV-l) samples prior to the data window and the M(N-l) samples 
20 after the data window. 

2) Apply the initial MMSE filter from (13) to obtain the initial impulse response 
estimates [x,{-(M -l)(AM)}, ••-,*,{- 1}, ^{0}, • • • , Jt,{L}, • ■ , x x {L-l-±{M -1)( AM)}] . 

3) Compute the initial power estimates A(0 = |*i(O| 2 for I = -(Af-1)(AM),-- S L-1-KAM)(AM) 
which are used to compute the filters w,(£) as in (11), and then apply to y(£) to 

25 obtain [i 2 {-(M-2)(AM)},.- s * 2 ^ 

4) Repeat 2) and 3), changing the indices where appropriate; until the desired length-L data window 
is reached. 

The initial estimate of the impulse response found by applying the MMSE filter is used as a priori 
information to reiterate the MMSE filter and improve performance. This is done by employing the MMSE 

30 filter formulation from (11) in which the respective powers of the impulse response coefficients are taken 
from the current estimate. It has been found that two or three reiteration steps allow the RMMSE filter to 
exclude the effects of scatterers prior to the data window, as well as to suppress the spatial ambiguities very 
close to the level of the noise floor. The RMMSE filter does especially well when the impulse response is 
somewhat sparsely parameterized (i.e. highly spiky), as is the case for system designed to have high 

35 resolution. Note that each reiteration step will reduce the number of coefficient estimates by 2(N -l). To 
counteract this, it is necessary to increase the size of the data window by 2M(N -l) samples, where M is 
the number of reiteration steps. Typically, however, L» N so that this reduction in data window size is 
negligible. 

An important factor regarding practical implementation is the non-singularity of the Afx/V matrix 
40 ([C(0]+ [R]) • This can be addressed by instituting a nominal level for which the estimated coefficients are 
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5 not allowed to go below. An alternative to this would be to re-estimate only those coefficients that are above 
some threshold since the small-valued coefficients do not contain a detectable scatterer. 

To demonstrate the performance of the RMMSE algorithm we examine three cases. The first case is 
typical of the scenario often addressed for LS estimation techniques and consists of an impulse response 
with a single large scatterer in noise and low-level clutter. The second case we address is when there is a 
10 second large scatterer that resides just prior to the data window. In this case it is expected that the LS 
estimator will substantially degrade since this important region is not accounted for in the received signal 
model. The final case is a more realistic scenario in which there are several spatially dispersed scatterers of 
varying magnitude. 

For the first two cases, the waveform selected is the length N = 30 polyphase modulated Lewis- 
15 Kretschmer P3 code, which on receive (after down-conversion to baseband in the receiver front-end) is 
defined as 



s(n) = exp 



n = 0X-,N-l. (15) 



20 The noise and clutter are modeled as zero-mean Gaussian with powers set to 70 dB and 50 dB below the 
signal power, respectively. For both cases we perform two iterations of the RMMSE algorithm and then 
compare the RMMSE-esti mated impulse response with the true impulse response, as well as with the results 
obtained from using LS and the matched filter. The impulse response for the first case is depicted in Figure 
3 for a single scatterer present in the data window which consists of 200 range gates. Figure 4 illustrates the 

25 results from the different estimation techniques in which the matched filter experiences significant spatial 
ambiguities while both LS and RMMSE have suppressed the spatial ambiguities so that the true scatterer 
location is evident and smaller nearby scatterers could be detectable. 

For the scenario just described, the LS and RMMSE estimators perform almost identically. However, 
when there is a significant scatterer present just prior to the range window, as depicted in Figure 5, the LS 

30 estimator is expected to degrade substantially. From Figure 6, the LS estimator truly does suffer from 
severe mis-estimation when a significant scatterer cannot be expressed in the model. However, the 
performance of the RMMSE estimator is indistinguishable from the previous case where only a single 
scatterer was present. This is obviously due to the fact that the RMMSE estimator takes the previous 
estimates of the surrounding impulse response coefficients into account when estimating the impulse 

35 response. 

Table 1 presents the Mean Squared-Error (MSE) for both cases discussed using LS, RMMSE, and 
normalized matched filtering. The MSE is averaged for all 200 coefficients in the data window and over 100 
runs using the same scatterer(s), with the clutter and noise distributed according to a zero-mean Gaussian 
distribution for each run independently. 

40 
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Table 1. MSE performance comparison 





Norm. MF 


LS 


RMMSE 


Case 1 


-123 dB 


-40.9 dB 


-35.8 dB 


Case 2 


-11.2 dB 


-11.3 dB 


-34.7 dB 



For case 1 (single scatterer), the normalized matched filter has a large MSE due to spatial ambiguities, 
while LS and RMMSE perform nearly the same with LS just marginally better. However, for case 2 
10 (scatterer just prior to the data window), the MSE attained by LS degrades to nearly that of the normalized 
matched filter, while the RMMSE maintains roughly the same MSE as in the previous case. 

For the third case, we examine the performance for the three estimation techniques when the impulse 
response contains several dispersed scatterers of varying magnitude (both within and outside of the data 
window). Figure 7 illustrates the significant improvement in estimation performance of the RMMSE 
15 estimator as compared to matched filtering and LS estimation. Most notable is the resolvability of two very 
closely spaced scatterers. RMMSE completely resolves the two scatterers whereas for LS and matched 
filtering it is not completely clear whether it is one or two scatterers. 

Obviously many modifications and variations of the present invention are possible in the light of 
the above teachings. It is therefore to be understood that the scope of the invention should be determined by 
20 referring to the following appended claims. 
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